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ABSTRACT 

We present results from a new XMM-Newton observation of the high-redshift quasar 
RXJ1028. 6-0844 at a redshift of 4.276. The soft X-ray spectral flattening, as reported 
by a study with ASCA previously (Yuan et al. 2000, ApJ 545, 625), is confirmed to 
be present, however, with reduced column density when modelled by absorption. The 
inferred column density for absorption intrinsic to the quasar is 2.1(^Q'g) x 10 22 cm~ 2 
for cold matter, and higher for ionised gas. The spectral flattening shows remarkable 
similarity with those of two similar objects, GB 1428+42f 7 (Worsley et al. 2004, MN- 
RAS 350, L67) and PMN J0525-3343 (Worsley et al. 2004, MNRAS 350, 207). The 
results improve upon those obtained from a previous short-exposure observation for 
RXJ1028. 6-0844 with XMM-Newton (Grupe et al. 2004, AJ 127, 1). A comparative 
study of the two XMM-Newton observations reveals a change in the power-law photon 
index from T ~1.3 to 1.5 on timescales of about one year. A tentative excess emis- 
sion feature in the rest-frame 5-10 keV band is suggested, which is similar to that 
marginally suggested for GB 1428+4217. 
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1 INTRODUCTION 

High-redshift quasars are powerful cosmological probes to 
study the evolution of massive black holes and quasar en- 
vironments in the early universe. Previous X-ray observa- 
tions suggested the presence of soft X-ray spectral flatten- 
ing 1 in some radio-loud quasars at redshifts z—2~3 (Wilkes 
et al. 1992, Elvis et al. 1994, Cappi et al. 1997, Fiore et 
al. 1998, Yuan & Brinkmann 1998, Reeves & Turner 2000). 
This result is strengthened and extended to higher red- 
shifts by its detection in a few extremely X-ray /radio-loud 
quasars at z >4, namely RX J1028. 6-0844 (Yuan et al. 2000) , 
GB 1428+4217 (Boiler et al. 2000, Fabian et al. 2001b), 
and PMN J0525-3343 (Fabian et al. 2001a), with ROSAT, 
ASCA, and BeppoSAX. These objects seem to have charac- 
teristics typical of blazars (Fabian et al. 1997, 1998; Zickgraf 
et al. 1997, Moran & Helfand 1997, Hook & McMahon 1998). 
The most plausible explanation for this effect is photoelec- 
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trie absorption of soft X-rays by associated medium with 
column densities of 10 22_23 cm -2 , although intrinsic spec- 
tral flattening cannot be excluded. The physical implication 
of this effect has been discussed extensively in the literature 
in terms of excess absorption (Elvis et al. 1998, Yuan et al. 
2000, Fabian et al. 2001a, b) and intrinsic breaks in the X-ray 
spectra of blazars (Fabian et al. 2001a, b). 

The contemporary X-ray observatories XMM-Newton 
and Chandra should be able to test these previous findings. 
Indeed, soft X-ray spectral flattening has been confirmed to 
be present in GB 1428+4217 and PMN J0525-3343 (XMM- 
Newton, Worsley et al. 2004a, b), and in some other objects 
at lower redshifts, e.g. PKS 2126-0158 at z=3.27 {XMM- 
Newton, Ferrero & Brinkmann 2003; BeppoSAX, Fiore et al. 
2003). Tentative evidence was also found in the combined 
spectra of several z > 4, moderately radio-loud quasars 
(Chandra , Bassett et al. 2004). 

RX J1028. 6-0844 was first detected as an X-ray source 
in the ROSAT All-sky Survey (RASS) and was identified as 
a quasar at z=4.276 (Zickgraf et al. 1997). It is also a ra- 
dio source (PKS B 1026-084) with a flux density of 220 mjy 
at 5 GHz (Otrupcek & Wright 1991) and a flat radio spec- 
trum. Its X-ray colors in the ROSAT energy band imply a 
hard spectrum (Zickgraf et al. 1997). Its first X-ray spec- 
trum, as obtained with a long ASCA observation made 
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in 1999, flattens substantially towards soft X-ray energies 
(Yuan et al. 2000); an excess (cold) absorption model re- 
quired a column density of ~ 2 x 10 22 cm -2 for local ab- 
sorber or ~ 2 x 10 23 cm -2 for absorber intrinsic to the 
quasar. A later short-exposure XMM-Newton observation 
found, however, only marginal evidence for excess absorp- 
tion (Grupe et al. 2004). In this paper we report on a new 
XMM-Newton observation of RX J1028.6-0844 with an ex- 
posure much longer than the previous observation. The mea- 
sured X-ray spectrum — with substantially improved photon 
statistics — confirms the presence of the soft X-ray spectral 
flattening as detected by ASCA. The observations and data 
reduction are described in Sect.[5] We present the spectral 
analysis in Sect.|21 including a re-analysis of the previous 
XMM-Newton observation. Discussion of the results is given 
in Sect.[I] including comparisons with previous observations 
and with other similar objects. Conclusions are summarised 
in Sect. GO We adopt H =71 kms" 1 Mpc -1 , r* A =0.73, and 
f2 m =0.27. The Galactic column density in the direction of 
RXJ1028.6-0844 is N§ al = 4.59 x 10 20 cm -2 (Dickey & 
Lockman 1990). Errors are quoted at the 1 a level for one 
parameter of interest unless stated otherwise. 



2 OBSERVATION AND DATA REDUCTION 

The quasar RXJ1028. 6-0844 was observed with XMM- 
Newton on June 13th, 2003 during satellite revolution 643 
(observation ID 0153290101). The EPIC (European Photon 
Imaging Camera) MOS1, MOS2, and PN cameras were op- 
erated in the 'primary full window' imaging mode and a thin 
filter was used to screen out optical/UV light. The observa- 
tional log is shown in Table 1. The XMM-Newton Science 
Analysis System (SAS, v. 6.0) and the most up-to-date cali- 
brations (August 2004) were used for data reduction. We fol- 
lowed standard data reduction and screening procedures. A 
fraction of the observation period suffered from high flaring 
background caused by soft protons. By inspecting the light 
curve of energy E >10keV, single events in the whole field 
of view, these periods were identified as having count rates 
higher than lctss" 1 and 0.35 ctss -1 for PN and MOS de- 
tectors, respectively, as recommended by the XMM-Newton 
Science Operation Centre (SOC). 

The quasar was detected at a sky position 
RA=10 h 28 m 38 a .84, Dec=-08° 44' 38". 3 (J2000), 0".6 
away from the position of its radio counterpart (Simbad 
database) and 2".l of its optical counterpart (Zickgraf 
et al. 1997). Source X-ray events were extracted from a 
circle of 32 "radius, which corresponds to the ~ 87 per 
cent encircled energy radius. Background events were ex- 
tracted from source-free regions using a concentric annulus 
of 52/128 "radii for the MOS detectors, and circles of 
32 "radius at the same CCD read-out column as the source 
position for the PN detector. X-ray images, light-curves, 
and spectra were generated from the extracted, cleaned 
events for the source and background. Parameters for data 
screening and source extraction are listed in Table 1 . No 
photon 'pile-up' problem was found, as expected, nor effect 
of low-energy noise above 0.2 keV. We used the spectral 
range of 0.2-10 keV for PN and 0.3-10 keV for MOS. The 
EPIC response files (rmf and arf) were generated using 



the source information on the detectors. The EPIC spectra 
were re-binned to have a minimum of 30 counts in each bin. 

The source profile in the 0.3-2 keV band was compared 
against the point spread function of the detectors (FWHM 
of 5 "for MOS and 6 "for PN) and was found to be con- 
sistent with a point-like source. There was no significant 
variability found during the 40 ks exposure, though a ~10 
per cent drop in count rates (averaged over ~5ks) was 
marginally detected, from 0.316±0.09 cts s _1 at the begin- 
ning to 0.271±0.08 cts s _1 towards the end of the observa- 
tion. 



3 X-RAY SPECTRAL ANALYSIS 

Due to an increase in the surface charge loss properties of 
the CCDs, which degrades the energy resolution, there has 
been a time-dependent, significant change in the low-energy 
redistribution properties of the MOS cameras. This effect 
has been taken into account in the most up-to-date MOS 
calibration files; however, some small systematic uncertainty 
may still remain at below 0.5 keV (Kirsch 2004, Kirsch et al. 
2004). To minimise possible biases induced in the results, 
we treat the MOS spectra in two ways and compare the 
results. The first was to simply omit the spectral range below 
0.5 keV; the second was to use the spectral range down to 
0.3 keV and introduce a systematic error of 2 per cent in the 
0.3-0.5 keV band (as recommended in Kirsch 2004, Kirsch 
et al. 2004). As seen in Tabled these two methods yielded 
statistically consistent results. We thus formally quote the 
results obtained using the 0.3-10 keV band. XSPEC (v.11.3) 
was used for spectral fitting. 



3.1 Soft X-ray spectral flattening 

3.1.1 Power-law model with local cold absorption 

We began by fitting the spectra of each detector individually 
with a single power-law model 2 modified by neutral absorp- 
tion with a column density Nh as a free parameter. This 
model gave acceptable 3 fits to the PN and the MOS1 spec- 
tra (see Table|5J, but not to the MOS2 spectrum [reduced 
X 2 =1.4 for 53 degree of freedom (d.o.f.), i.e. a null hypothesis 
probability of P nu n=0.02 only]. Inspection of the \ 2 residu- 
als of the fit singled out one energy bin at 1.82 keV (a width 
of 80 eV) , which contributed 1 1 out of the total \ 2 of 76. The 
feature can be fitted with either a narrow notch feature at 
£=1.83±0.02 keV (a width of 49±?| eV and a covering frac- 
tion O.99I0 53) or a Gaussian absorption line (£=1.83 keV) 
of infinitely small width. This energy corresponds to 9.7 keV 
in the quasar's rest-frame, at which no known physical ab- 
sorption feature is present. On the other hand, it is coinci- 
dent with the instrumental absorption feature at 1.84 keV 
due to the Silicon edge. Furthermore, it appears in neither 
the PN nor the MOS1 spectrum. Excluding the 1.82 keV 
bin reduced the x 2 by 12, and resulted in an acceptable fit 



2 E p ^ Q {E) °c E~ r , where Y is the photon index. 

3 We regard a model fit as acceptable if the null hypothesis proba- 
bility derived from the fit is greater than 10 per cent (P nu n >0.1), 
as is commonly quoted. 
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Table 1. Summary of the XMM observations and data reduction information for each EPIC detector 





PN 


MOS1 


i vi v_y o — 


observation duration (ks) 


41.5 


43.1 


43.1 


good exposure (ks) 


15.8 


21.3 


22.1 


energy band used (kcV) 


0.2-10 


0.3-10 


0.3-10 


events pattern used 


0-4 


0-12 


0-12 


source extraction radius 


32" 


32" 


32" 


source+BGD counts 


5024 


2006 


2080 


net source counts 


4854 


1920 


1986 


source count rate (10 — 2 ctss _1 ) 


30.7 ± 0.5 


9.0±0.2 


9.0±0.2 


fiux a (lO-^ergs -1 cm -2 ) 


1.16±0.04 


1.14±0.07 


1.14±0.06 


flux 6 (lCr^ergs" 1 cm" 2 ) 


1.17±0.04 


1.15±0.07 


1.15±0.06 



a Flux in the 1-10 keV band; the model of a power-law with intrinsic absorption at 2=4.276 plus 7V^* al is used. 
b Galactic absorption corrected flux in the 1—10 keV band; the same model as in a is used. 



(x 2 =l-2 for 52 d.o.f.) with fitted parameters in good agree- 
ment with those for MOS1 within 1-a errors. We therefore 
consider the disagreement between MOS2 and MOS1/PN to 
be due to inadequate MOS2 calibration around the Silicon 
edge, and ignore this energy bin hereafter. Parameters (in- 
cluding normalisations) for the two MOS spectra were tied 
together to be the same in joint fitting (MOS1+2). The fit- 
ted Nh of ~ 1.1 x 10 21 cm -2 ) is significantly higher than 
iVH al (TableEJ. The photon index is now F ~1.55, typical of 
blazars. 

A power-law with fixed Galactic absorption 
(4.59 x 10 20 cm -2 ) yielded unacceptable fit and a flat 
photon index (r=1.3-1.4, Table|5J. The improvement in 
X 2 for the fit with freely fitted Nh over that with fixed 
iVH=-/Vjf al is substantial — \ 2 was reduced by 36 for PN 
and by 31 for a joint MOS1+2 fit for one additional free 
parameter. Applying the .F-test 4 (Bevington & Robinson 
1992) gives a chance probability <C0.001. 

For fixed Nh=N§^ 1 , acceptable fit was obtained only 
within the restricted 1-10 keV range. We plotted in Fig.Q 
the data, the best-fit model (for a joint PN and MOS 1+2 
spectral fit) and its extrapolation down to the low-energy 
end of the detectors, and the data-to-model ratio as resid- 
uals. A systematic deficit of photons below 1 keV is clearly 
indicated. Fig.|5| shows the confidence contours for the free- 
fitted total Nh (Galactic plus excess) and F for PN and 
MOS1+2, respectively. Absorption in excess of N§^ is ev- 
ident. It is noted that PN gives systematically lower Nh 
values than MOS. We regard the results from MOS to be 
more reliable than those from PN in consideration of the 
EPIC cross calibration uncertainties as discussed in details 
in AppendbJXl 



4 Recently, Protassov et al. (2002) have questioned the validity of 
using the F-test to test for the significance of adding an additional 
spectral component, as this involves testing a hypothesis that is 
on the boundary of the parameter space. We note that the .F-test 
we used here is not affected by the boundary condition problem, 
because we are testing for the significance of the different values 
in JVjj between Nh=N^ 31 (null hypothesis) and the free-fitted 
value, rather than for the addition of a model component. 
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Figure 1. The spectra of the PN (circles), MOS1 (squares), and 
MOS2 (stars) cameras and the residuals as data-to-model ratio. 
The model is the best-fit power-law (with Galactic absorption) 
to the joint PN+MOS1+MOS2 spectra within the restricted 1- 
10 keV energy band and is extrapolated down to the low energies. 
A systematic deviation from a power-law model with Galactic 
absorption towards low energies (below 1 keV) is clearly indicated. 

3.1.2 Absorption intrinsic to the quasar 

This is the most plausible postulate in consideration of 
the statistical argument presented in previous work (see 
Introduction for references). For neutral absorber, Nh — 
2 x 10 22 cm" 2 was found assuming cosmic abundances. The 
confidence contours for the excess Nh versus V are shown 
in Fig.dl for PN (dashed) and MOS1+2, respectively. If 
the metallicity of the absorber is lower than the cosmic 
value, which is not unexpected at such a high redshift — 
could be about 10 per cent or less (e.g. Lu et al. 1996, Pet- 
tini et al. 1997, Prochaska & Wolfe 2000), the iV H would 
be correspondingly higher than the value give here. The 
absorption-corrected luminosity in the quasar rest-frame is 
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Table 2. Results of X-ray spectral fits 



Detector 



X 2 /dof 



pb 

■'null 



local neutral absorption, free Nn 
PN C 

MOS1 (0.5-10keV) 
MOS1 (0.3-10keV) 
MOS2 (0.5-10keV) 
MOS2 (0.3-10keV) 
MOS1+2 (0.5-10keV) 
MOS1+2 (0.3-10keV) 



8.7±0.7 


1.53±0.03 


137.8/147 


0.69 


12.8±2.5 


1.54±0.06 


57.7/52 


0.27 


12.7±2.0 


1 54+O.O6 


68/67 


0.44 


11.3±2.5 


1 C-+0.06 
1,0 ' -0.03 


63.7/52 


0.13 


9.0±1.6 


1.53±0.05 


76/67 


0.21 


12.1±1.7 


1 5fi+ - 05 


122.8/107 


0.14 


11.1±1.2 


1.54±0.04 


150/137 


0.21 



local neutral absorption, 


fixed Nn 










PN 


4.59 


(fix) 


1.37±0.02 


173/148 


0.07 


MOS1 (0.5-10keV) 


4.59 


(fix) 


1.36±0.04 


69/53 


0.07 


MOS1 (0.3-10keV) 


4.59 


(fix) 


1.30±0.03 


91/68 


0.03 


MOS2 (0.5-10keV) 


4.59 


(fix) 


1.43±0.04 


71/53 


0.05 


MOS2 (0.3-10keV) 


4.59 


(fix) 


1.40±0.03 


84/68 


0.09 


MOS1+2 (0.5-10keV) 


4.59 


(fix) 


1.40±0.03 


141/108 


0.017 


MOS1+2 (0.3-10keV) 


4.59 


(fix) 


1.35±0.02 


181/138 


0.008 



power-law + 7V^* al and excess 


absorption 


at z = 4.276 






Detector 


jyexc 


r 


X 2 /dof 


-Pnull 


PN 

MOS1+MOS2 (0.5-10keV) 
MOS1+MOS2 (0.3-10keV) 


149±30 
393+99 

222±46 


1.49 ± 0.03 

i 53+O.O4 

1.49±0.03 


136/147 
122/107 
153/137 


0.15 
0.17 



broken power-law + Galactic absorption (PN+MOS1+2) 







break 


Tlow-E 


X 2 /dof 


-Pnull 


PN 




o 5n+ - 04 


-1.55+ - 70 


133/146 


0.77 


MOS1- 


k2 (0.5-10keV) 


1 1+ 01 


0.86±°;» 


122/105 


0.12 


MOS1- 


•-2 (0.3-10keV) 


1 1+0-2 

1 - L — o.i 


0-92th 


148/135 


0.21 



a Column density of hydrogen in units of 10 20 cm -2 
b Null hypothesis probability of the x 2 test for the fit 

c The A^h values from PN are possibly under-estimated; see Sect, l4,l"Tl and AppendixlAl for discussion. 

d The break energy of broken power-law in keV. When no errors are given, the value is unconstrained within a physically meaningful 



9.2 xlO^ergs" 1 in 1-10 keV and 2.64 x 10 47 erg s" 1 in 1- 
50keV, respectively. 

If the absorber is close enough to the central source, 
the gas is likely to be ionised. Indeed, the optical spectra of 
RX J1028.6-0844 taken by Peroux et al. (2001) and Zickgraf 
et al. (1997) show no significant Lyman limit absorption at 
912 A . An estimate of the optical depth of the Lyman limit 
absorption places an upper limit on the neutral hydrogen 
column density to be Nn i < 10~ 17 cm -2 along the line-of- 
sight to the optical-UV emission region (Yuan et al. 2000). 
If the X-ray absorber also intercepts the optical-UV light, at 
least moderate ionisation of the gas is required. The lack of 
strong optical-UV extinction can be explained by ionised, 
dust-free absorber. It is interesting to note that such an 
optical-UV property seems to be common among this type 
of objects (Yuan et al. 2000, Fabian et al. 2001a,b). We tried 
to model the excess absorption with ionised absorption (ab- 
sori in XSPECT). The ionisation parameter (as defined in 
Done et al. 1992), however, could not be constrained with 
the current data ranging from an almost neutral to highly 
ionised absorber. The results of such an analysis are given in 
Fig.[3| where the data from the two XMM-Newton observa- 



tions are combined in order to achieve a better constraint on 
the parameters (see Sect . 13. 2^21 . With increasing ionisation 
state, the required column density increases from ~ 2 x 10 22 
to ~ 1 X 10 23 cm" 2 . 

No redshifted K-shell absorption edges from iron ions 
(Scdgc=7. 1-9.3 keV in the rest-frame from Fei to Fexxvi) are 
detected. Assuming cosmic abundance of iron (3.4 x 10 -5 ), 
a K-shell edge with optical depth r ~ 0.02-0.1 is pre- 
dicted from the above Nn range for neutral-highly ionised 
iron. This is consistent with the upper limit on r esti- 
mated from the joint PN and MOS1+2 spectra, which 
ranges from 0.1 to 0.2 (90 per cent level) for ions from 
Fei to Fexxvi (absorption cross section from 3.8 x 10 -20 
to 3.3 x 10" 20 cm 2 atom" 1 ). 



3.1.3 Intrinsic spectral break 

We also considered the possibility that the soft X-ray spec- 
tral flattening is an intrinsic feature. The physical implica- 
tion of a break in the intrinsic X-ray spectra of blazars has 
been discussed in Fabian et al. (2001a,b) in the context of a 
cut-off in the energy distribution of electrons or of soft seed 
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N H (10 22 cm" 2 ) 



Figure 2. Confidence contours of fitted total column density 
and photon index for the model of a power-law with local neu- 
tral absorption. The contours are at the 68, 90, and 99 per cent 
confidence, respectively, for two interesting parameters. Solid con- 
tours: the joint MOS spectra in 0.3— lOkeV; dashed contours: the 
PN spectrum. Also indicated by lines are the Galactic column 
density (dashed) and its conservative 30% uncertainty range (dot- 
ted). Absorption in excess of the Galactic value is evident for such 
a model. 
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Figure 3. Column density of intrinsic absorber versus power-law 
photon index as measured with MOS (solid) and PN (dashed) 
detectors. The absorber is assumed to be at the quasar redshift 
4.276 and to have cosmic abundances. The contours correspond to 
confidence levels of 68, 90, and 99 per cent. The results from the 
current (2003) and the previous XMM observations (2002) are 
plotted. The variation in the spectral slope is significant, while 
there is no change in the absorption column density. 

photons for Compton scattering. We modelled the spectral 
flattening with, as an approximation, a broken power-law 
modified by local absorption. Acceptable fits could be ob- 
tained for both PN and MOS1+2. However, Ah and the 
low energy photon index ri ow _E could not be constrained 
due to strong coupling. We thus fixed N H =N§ al . The fit- 
ted high-energy photon indices Fhigh-E are 1.44±0.04 (PN) 
and 1.49±0.06 (MOS1+2), and the break energy E hlcak and 
low-energy index Piow-E are listed in Tabled This model 
gave acceptable fits which are statistically indistinguishable 
from the models of power-law with either local or intrinsic 
absorption. 



3.2 Comparison with a previous observation and 
spectral variability 

3.2.1 A previous XMM-Newton observation 

RXJ1028.6-0844 was previously observed with XMM- 
Newton in revolution 445 with a short exposure of 7ks in 
May, 2002 (PL S. Mathur, observation ID: 0093160701). 
The results, as published in Grupe et al. (2004), gave ab- 
sorption A r H values similar to what we find here, but, a 
flatter spectral slope of T ~ 1.3. In order to achieve a self- 
consistent comparison of the two observations — free from 
the effects introduced by different versions of the evolving 
calibration and data processing software, we also analysed 
the data from that observation. The data were taken from 
the XMM-Newton science archive. The observation was de- 
scribed in Grupe et al. (2004). We used exactly the same 
data screening criteria and source/background extraction 
regions as used for the current observation (see Sect and 
TableQ. The good exposure and source count rate are 3.7 ks 
and 0.39 ctss" 1 for PN (0.2-10 keV), and 6.9/7.0 ks and 
0.12/0.13ctss" 1 for MOS1/2 (0.3-10 keV). A comparison 
with Tabled reveals that the broad-band count rates were 
higher in the 2002 observation than in 2003 by about 30 per 
cent. 

The PN and MOS spectra were binned to have 
a minimum of 25 and 20 counts in each bin, respec- 
tively. The results of the spectral fits are in good agree- 
ment with those obtained by Grupe et al. (2004). For 
an absorbed power-law model, the total local absorption 
A H is 10.1(±i;°) xl0 20 cm- 2 for joint MOS1 and MOS2 
(MOS1+2) spectra and 6.0(tJ; ) x 10 20 cm -2 for PN. Again, 
the fitted Njj. is systematically lower for PN than for MOS, 
as found in the 2003 observation. While the absorption Ah is 
in good agreement between the two observations, the pho- 
ton indices are not. The photon indices obtained for the 
2002 observation are r 2 oo2=l-30 ± 0.06 and 1.27 ± 0.05 for 
MOS1+2 and PN, respectively, i.e. the spectrum was steeper 
during the observation of 2003 (r 20 03 = 1-53 ± 0.03). This 
result can be seen clearly in Fig-El where excess (neutral) 
absorption Ah is plotted versus V for the 2002 observation, 
assuming absorption is intrinsic to the quasar at z=4.276. 
The spectral steepening remains even when only the hard 
band 2-10 keV spectra were considered (r 2 oo2=1.33t ;?f 
and 1.23t u ;?3 for MOS1+2 and PN, respectively). The 
Galactic absorption corrected flux in the 1-10 keV band is 
1.9 xl0~ 12 ergs" 1 cm -2 (averaged MOS1 and MOS2 value). 

3.2.2 Joint spectral fit of the two observations 

We quantified the spectral variability by fitting jointly the 
spectra of the two observations. We used the MOS spectra 
only, in consideration of possible PN calibration uncertain- 
ties (see AppendixlAl . For MOS1 and MOS2 spectra from 
the same observation, all parameters (including normalisa- 
tion) were tied together. The results are summarised in Ta- 
bleU 

Firstly, since the fitted Ah of the two observations are 
in good agreement, we assumed that there was no variability 
in the absorption. Thus the Ah values for the two observa- 
tions were tied together in the fitting. Absorbed power-law 
models were used. As a test, we tied r2oo2=F2oo3 together in 
the fitting, which resulted in a fit only marginally acceptable 
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(x 2 =242 for 219 d.o.f.). Setting the two indices as indepen- 
dent parameters improved the fit significantly, reducing \ 2 
by 24 for one additional free parameter; the F-test [Beving- 
ton & Robinson 1992; see Footnote^] for the argument for 
the validity of the F-test used here regarding the bound- 
ary condition warning discussed by Protassov et al. (2002) .] 
gives a chance probability <C 0.001 for F2oo2 and F2oo3 being 
the same. The fit is good (x 2 =217 for 218 d.o.f.), indicating 
that the spectra of the two observations can be fitted well 
with two different continua (in slope and normalisation) at- 
tenuated by the same amount of absorption. The fluxes in 
the 0.2-1 keV band were comparable in the two observa- 
tions, while in the 1-10 keV band it was higher by a factor 
of two in May 2002 compared to June 2003 (TableOJ. The 
two power-law continua, before attenuation by any absorp- 
tion, cross over at ~ 0.4 keV in the observer's frame, i.e. 
~ 2 keV in the quasar rest-frame. 

The combined data set improves the spectral photon 
statistics and gives a better constraint on the excess ab- 
sorption. For neutral absorber intrinsic to the quasar, Ah 



=2.1( 



+° :t)xio' 



cm , i.e. excess absorption is evident. An 



ionised absorber model also gives a good fit, suggesting that 
the ionisation status is unconstrained. The confidence con- 
tours for the ionisation parameter and Nb. are shown in 

Fig.m 

Secondly, we tested whether the variation in spectral 
slope could result from variability in the absorption, such 
as from the ionisation parameter. We fitted the ionised ab- 
sorber model with independent ionisation parameters and 
Ah, but tied the indices together (F2oo3=r2oo2). No satis- 
factory result could be obtained for this model; the best fit 
was significantly worse than that with the above variable-F 
model (Ax 2 = 14). This is not surprising, as the difference 
in r arises primarily in the hard 1-10 keV band, which cor- 
responds to 5.3-53 keV in the quasar rest-frame. At such 
high energies the amount of absorption of X-ray photons is 
decreasing dramatically. 

We also fitted a broken power-law to the data with ab- 
sorption fixed at A r ^f al . There is a marginal indication of a 
higher break energy 1.6±0.2keV for the flatter spectrum in 
May 2002 than for the steeper spectrum in June 2003; how- 
ever, the significance is low. The low-energy photon indices 
and the normalisations at 1 keV (below or close to the break 
energies) are comparable in the 2002 and 2003 observations. 

No spectral variability is detected within a 40 ks dura- 
tion in the June 2003 observation with XMM-Newton. 



4 DISCUSSION 

4.1 The presence of soft X-ray spectral flattening 

4.1.1 XMM-Newton results 

We have shown the presence of soft X-ray spectral flattening 
in the z=4.276 quasar RX J1028. 6-0844 using an observation 
made with XMM-Newton. The result confirms the previous 
report based on ASCA data (Yuan et al. 2000). In the ex- 
cess absorption scenario, the derived absorption Nh from 
our analysis are consistent with those obtained by Grupe et 
al. (2004) from a previous short XMM-Newton observation. 
In that study the authors suggested that strong excess ab- 
sorption was marginal. This is not surprising given the lower 



Table 3. Joint fits to the current (2003) and a previous (2002) 
XMM-Newton observations. Only the MOS spectra were used. 
The absorption 7V"h is assumed to be the same in the two obser- 
vations and the photon index is freely fitted. 



parameter 



June 2003 



May 2002 



power-law with neutral absorption 
photon index 

normalisation @lkeV (10 — 4 ) 
total absorption iVjj (10 22 ) 
flux" 0.2-1/1-10 keV (10 -13 ) 

neutral absorption at z=4-276, fixed 
photon index 

normalisation @lkeV (10 — 4 ) 
excess absorption TVjj (10 22 ) 
flux 6 0.2-1/1-10 keV (10~ 13 ) 
luminosity 1-10/1-50 keV (10 47 ) 

broken power-law, fixed 7V^ al x' 2 
break energy (keV) 
low-energy index 
high-energy index 
normalisation @lkeV (10 -4 ) 



1.53±0.04 1.31±0.04 
1.77±0.07 2.13±0.09 

O.llztO.Ol (tied) 
1.2/11.5 1.4/19.1 



,+0.4 
-0.3 



1.48±0.03 

1.65±0.05 
2.F 
1.7/11.6 
0.92/2.58 

/d.o.f.=214/215 

1.1±0.1 1.6±0.2 

0.92±0.11 0.87±0.09 

1.50±0.04 1.33±0.07 

1.54±0.06 1.75±0.06 



X 2 /d.o.f.=221/218 
1.26±0.04 
1.99±0.07 
(tied) 
1.8/19.2 
1.06/3.88 



a Galactic absorption uncorrected flux in units of ergs -1 cm -2 . 
b Galactic absorption corrected flux in units of ergs -1 cm -2 . 
c Absorption corrected luminosity in the quasar rest-frame in units 
of ergs . 




Figure 4. Confidence contours (at the 68, 90, and 99 per cent 
level) for the ionisation parameter and column density of the 
ionised absorber model. The result is obtained from a joint spec- 
tral fit to the MOS spectra from the two XMM-Newton observa- 
tions in 2003 and 2002 (see Sect.| 



signal-to-noise of their data and consequently the weaker 
constraints on Nh compared to this work, which benefits 
from a much longer exposure. 

We note that PN tends to give systematically lower Ah 
values than MOS. This is most likely due to discrepancies 
in the calibration below 1 keV between PN and MOS, as re- 
ported in the most recent XMM calibration status (Kirsh et 
al. 2004, see also XMM-Newton SOC XMM-SOC-CAL-TN- 
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0018 5 ). Although a complete solution is yet to be reached, 
preliminary indications suggest that PN is most likely the 
cause of the problem. A detailed discussion on this issue 
and a quantitative PN-MOS comparison taking the most 
up-to-date calibration into account is given in AppendbdAl 
In summary, we consider the Nn values derived from MOS 
spectra to be more reliable. Furthermore, it may also be 
the case that both PN and MOS yield systematically higher 
fluxes at low energies compared to XMM-Newton RGS and 
Chandra. If this turns out to be the case, the true column 
in the absorption scenario could be even higher than what 
is reported here. 

4-1-2 Comparison with previous results 

The total Ah assuming local absorption is ~ (0.11 ± 
0.01) x 10 22 cm -2 measured from the XMM-Newton obser- 
vations (joint MOS from the two observations). This value 
is a factor of 2-3 times smaller than the value measured 
from ASCA (Yuan et al. 2000). For intrinsic absorption at 
z = 4.286, the N-r inferred by XMM-Newton (a few times 
10 22 cm -2 ) becomes about 10 times smaller than that ob- 
tained by ASCA. Intrinsic variability in the absorption can- 
not be ruled out. However, we speculate that a systematic 
difference in the instrumental calibration of the two missions 
might play at least a partial role. This is because a similar 
trend was also found for GB 1428+4217 and PMN J0525- 
3343 (Worsley et al. 2004a, b). It is not clear which instru- 
ment causes the difference. In the cases of GB 1428+4217 
and PMN J0525-3343, the previous results before XMM- 
Newton were obtained by a joint fit of both the ASCA and 
BeppoSAX spectra (Fabian 2001a,b). This fact, together 
with the aforementioned XMM EPIC calibration issue, sug- 
gests that perhaps both missions, rather than merely XMM- 
Newton, may be the cause. Improved XMM EPIC calibra- 
tion and independent investigations by other instruments 
( Chandra or XMM-Newton RGS) are needed to resolve this 
problem. 

Comparing the observations of XMM-Newton in 2003 
and of ASCA in 1999, the spectral photon indices are con- 
sistent within their mutual 1-a errors; no significant flux 
variability (>10 per cent) is detected in the 1-10 keV band. 

4-1.3 Comparison with other objects 

We compare the spectral shape of RXJ1028. 6-0844 with 
those of GB 1428+4217 and PMN J0525-3343 in the quasar's 
rest-frame. Following Worsley et al. (2004b), we produced 
the data-to-model ratio for RX J1028. 6-0844 where the 
model is the best-fit power-law with Galactic absorption in 
the restricted 1-10 keV band. The result is plotted in Fig.|S| 
together with those for GB 1428+4217 and PMN J0525-3343 
from Worsley et al. (2004b, their Figs. 2 and 5). It should be 
noted that the data-to-model ratio is free from the effects 
of instrument response, Galactic absorption, and redshift. 
It can be seen that the three spectra agree remarkably well 
in terms of the break energy and the shape of the spectral 
cutoff. The Ah values of intrinsic (cold) absorbers of a few 
times 10 22 cm _2 as measured by XMM-Newton in these three 

5 http:/ /xmm. vilspa.esa.es/es/external/xmm_sw_cal/calib/indox. 
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Figure 5. Data-to-model ratio in the rest-frame for RX J1028.6- 
0844 (stars), where the model is the best-fit power-law with 
Galactic absorption in the restricted 1— lOkeV spectral range. 
Only the PN spectrum (May 2003 observation) is plotted, 
which is rebinned for demonstration. Also plotted are those for 
GB 1428+4217 (filled squares) and PMN J0525-3343 (open cir- 
cles) for comparison, which were taken from Fig. 5 in Worsley et 
al. (2004b). The spectra are plotted in their respective rest-frame 
energies. 



objects are in good agreement. The striking spectral simi- 
larity shared by these objects at different redshifts argues 
for a real soft X-ray flattening against instrumental effects, 
and suggests a common nature to this phenomenon. 

Another similarity lies in their optical-UV properties, 
which argue for a highly ionised, dust-free absorber model 
(Yuan et al. 2000; Fabian et al. 2001a,b; Worsley et al. 
2004a,b). It is interesting to note that excess absorption in 
several z > 4, moderately radio-loud quasars, as tentatively 
suggested by their combined Chandra spectra, has also simi- 
lar Nn values of a few times 10 22 cm -2 (Bassett et al. 2004). 



4.2 X-ray spectral variability 

It is worth stressing that while spectral variability is indeed 
typical of flat-spectrum quasars, the extremely flat spectrum 
in the 2-10 keV band during the 2002 observation is consis- 
tent only within the 2 a uncertainty range with the limiting 
value of F ~ 1.5 for a relativistic distribution of particles 
emitting via synchrotron and inverse Compton (in the sim- 
plest hypothesis) . If the production of such fiat high-energy 
spectra are confirmed, revision of the widely accepted emis- 
sion scenarios would be required. 

Although flux variability over short timescales is a dis- 
tinctive characteristic of blazar emission, no significant vari- 
ations have been detected within a single observation. It 
should be noted, however, that due to the high redshift of 
the source, the intrinsic timescale sampled by the observa- 
tion (~ 2 hr) might be too short to detect significant varia- 
tions (for powerful quasars doubling timescales of the order 
~several hour to a day might be more typical, e.g. 3C279, 
shtmlWherle et al. 1998). 
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4.3 Excess emission around 5— lOkeV? 

Worsley et al. (2004b) pointed out possible excess emission 
at energies around 5-10 keV in the quasar rest-frame spec- 
trum for GB 1428+4217. The evidence is only marginal. 
Interestingly, the same spectral structure also appears in 
RXJ1028.6-0844, as can be seen in Fig.El The similarity 
of the energy position of this feature in two objects at dif- 
ferent redshifts is remarkable. If this feature is real, it may 
come from an additional spectral component which is peaked 
around 5-10 keV in the rest-frame. We modelled the XMM- 
Newton spectra taken in May 2003 by adding a steep power- 
law component in the above spectral models, following Wors- 
ley et al. (2004b). Both the PN and MOS1+2 spectra were 
fitted jointly to improve the statistics. A steep photon in- 
dex is yielded for the second power-law, as F=2.3-4.4 for 
a model with local absorption and F=l. 9-7.2 for intrinsic 
absorption at z=4.276 (90 per cent confidence range for 
1 interesting parameter). These values are consistent with 
that obtained for GB 1428+4217, V ~1.8-2.6 (Worsley et 
al. 2004b). The model does improve the fit in the 5-10 keV 
(rest-frame) band, though the statistical significance is not 
high (Ax 2 — —7 for 3 additional free parameters). 

If this excess emission feature proves to be real, it might 
be the first evidence for the presence of emission originating 
from bulk comptonisation on the soft photon field through 
which the relativistic jet propagates (Begelman & Sikora 
1987). The detection of such a feature could carry key clues 
to the amount of (cold) leptons flowing in the jet (Sikora 
& Madejski 2000). The lack of its detection in the majority 
of objects so far remains a puzzle. However, in most cases 
the X-ray emission might be dominated by the non-thermal 
emission from relativistic particles, and thus the possibility 
of detecting such a component could be limited to cases of 
particularly high jet Lorentz factors (which would shift its 
peak up to high energies) and/or low-states/steep power-law 
of the non-thermal relativistic component (the latter case 
could of course be tested, in principle). Observations with 
even higher signal-to-noise than the present ones or stacking 
spectra from different sources might be a way to clarify the 
issue. 



5 CONCLUSIONS 

We have presented a new X-ray spectroscopic study of the 
high redshift (z=4.276) quasar RX J1028.6-0844 with XMM- 
Newton. The high signal-to-noise spectrum confirms the 
presence of the soft X-ray spectral flattening, which was re- 
ported previously with ASCA data (Yuan et al. 2000). This 
spectral feature can be modelled by either excess absorp- 
tion of the quasar X-rays or an intrinsic break at ~lkeV 
in the X-ray spectra of the source. In the absorption sce- 
nario, the derived column density for cold absorber intrinsic 
to the quasar is 2.1(^°'g) xlO 22 cm -2 . This value is compa- 
rable to those reported in two similar objects GB 1428+4217 
(2=4.72) and PMN J0525-3343 (2=4.4) from XMM-Newton 
observations. The remarkable similarity in the shape of the 
spectral cutoff among these objects at different redshifts ar- 
gues against instrumental effects as an origin, but rather 
argues for a common nature to the soft X-ray flattening in 
high-z blazars. In terms of the soft X-ray spectral flatten- 
ing, the results are consistent statistically with and improved 



upon those obtained from a previous short-exposure obser- 
vation for RXJ1028. 6-0844 with XMM-Newton (Grupe et 
al. 2004). A comparative study of the two XMM-Newton 
observations revealed a spectral steepening from V ~1.3 in 
2002 to r ~1.5 in 2003, and a consequent drop in flux in the 
hard energy band above ~lkeV. 

The derived columns from XMM-Newton observations, 
however, are reduced when compared with the previous 
ASCA results (Yuan et al. 2000). We speculate that this 
might be due to systematic instrumental effects, probably 
inherent in both missions. Future improved XMM-Newton 
EPIC calibration and independent investigations by other 
instruments (such as Chandra) are needed to resolve this 
issue. 

A tentative excess emission feature in the rest-frame 
5-10 keV band is suggested, which bares remarkable simi- 
larity to that marginally imprinted in the X-ray spectrum 
of GB 1428+4217 (Worsley et al. 2004b). 
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Table Al. absorbed power-law model fit to the PN spectrum 
using corrected effective area ARF calibration 



absorber 


N H a r 


X 2 /dof 




observation 


in June 2003 






local 


11.1±0.8 1.54±0.03 


137/147 


0.70 


local 


4.59 (fix) 1.29±0.02 


218/148 


1.410" 4 


z=4.276 


231±28 1.50 ±0.03 


134/147 


0.77 



a Column density of hydrogen in units of 10 20 cm 2 



10-15 per cent with respect to MOS and by 20-30 per cent 
with respect to RGS (Kirsh et al. 2004). It appears that 
the problem arises mainly from the PN calibration at low 
energies, although it is not clear whether the redistribution 
matrix or the effective area is deficient. This means that 
PN spectra tend to give less absorption column than MOS 
and RGS, a trend similar to what we find in this work (Ta- 
ble|21 and Table|3J , as well as those obtained by Grupe et 
al. (2004). We thus argue that the A^h values derived from 
MOS spectra are more reliable. 

If the effect is purely due to an effective area problem of 
PN, we can test the PN-MOS consistency by taking into ac- 
count the reported difference in the calibration. To simplify 
the treatment, we assumed that the PN effective area be- 
low 1 keV was under-estimated by a factor of f(E), which is 
energy-dependent. f(E) was estimated from Fig. 10 in Kirsh 
et al. (2004), in which the MOS/RGS normalisation factors 
with respect to PN is plotted versus energies. We adopted 
conservative values of /(0.3keV) ~1.4, /(0.6keV) ~1.15, 
and /(IkeV) ~1.0. f(E) at any other energies E within 
0.3-1.0 keV was interpolated using a binomial function. We 
corrected the PN effective area (ARF) by multiplying by 
f(E) at the corresponding energy. The PN spectral fitting 
was repeated with is corrected ARF and the results are listed 
in TableEU The fitted total N H is (11.1±0.8) 10 20 cm -2 , in 
excellent agreement with the MOS results, while V remains 
unchanged. We therefore conclude that the Ah fitted from 
the PN spectrum is likely to be under-estimated. More quan- 
titative and reliable estimation of A^h in PN spectra must 
await the completion of the PN and MOS calibration at the 
low energies (Kirsh et al. 2004). 



APPENDIX A: EFFECT OF CALIBRATION 
UNCERTAINTY ON THE RESULTS 

In the most updated reports on the XMM-Newton calibra- 
tion status, Kirsh et al. (2004, see also XMM-Newton SOC 
XMM-SOC-CAL-TN-0018 6 ) demonstrated significant dis- 
crepancies at low energies (below ~lkeV) between the EPIC 
PN and MOS cameras, and the EPIC and the RGS 7 , as well 
as overall differences in the XMM-Newton instruments and 
those of Chandra. While MOS and RGS agree with each 
other in general, PN gives higher fluxes below 0.7 keV by 



6 http: / / xmm.vilspa.esa.es /es/ external / xmm_sw_cal / calib / index. shtml 

7 The Reflection Grating Spectrometer on-board XMM-Newton. 



